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Abstract 

The hadronic vacuum polarization contribution to the anomalous magnetic 
moment of the muon is parametrized by using the quark-resonance model 
formulated in jlj. In this context a recent prediction obtained within the 
ENJL model Q can be affected by two additional contributions: the next to 
leading corrections in the inverse cutoff expansion and the gluonic corrections. 
Motivated by the necessity of reaching a highly accurate theoretical prediction 
of the hadronic contribution to the muon g-2, we study in detail both the 
effects. 

Short Title: The muon g-2 in the Quark-Resonance model. 
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Effective chiral Lagrangians a la Nambu-Jona Lasinio are a good theoretical 
framework to understand the hadronic vacuum polarization contribution to the 
anomalous magnetic moment of the muon, given by the diagram in Figure 1. 

Phenomenological estimates of a 1 ^ = (g h — 2)/2 are obtained from the best fit of 
the e + e~ — > hadrons total cross section a h {t) through the usual dispersion relation 



/ dt K(t) a h (t), (1) 
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whith the QED function K(t) given by: 

Kit) = f dx 2 f^~ X l, 2 - (2) 
Jo x l + (1 — x)t/m~ 

The most recent numerical estimates give the following values: 



• 7.07(.066)(.17) ■ 1(T 8 [3] 

• 6.84(.ll) • 1(T 8 m 

• 7.100(.105)(.49) • 1(T 8 [|, (3) 

where the first error is statistical and the second one is systematic. In what 
follows we give a theoretical picture to understand the above numerical values. 

is related to the renormalized hadronic photon self-energy 11^ (Q 2 ) through the 
following integral B B[ 
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n^(Q 2 ) is given in terms of the vector two point function Uy(Q 2 ) which has been 
extensively analyzed in refs. [I], [7|]: 

n^(Q 2 ) = E Qi (K(Q 2 ) - nir(o)) = \{rt v {Q 2 ) - n^o)), (5) 

i=u,d,s 

where Qi = (2/3, —1/3, —1/3) are the charges of the SU(3) flavour quarks u,d,s 
and the renormalized photon self-energy satisfies the constraint 11^(0) = 0. 

Because the typical momenta of the off-shell photons are of the order of the 
squared muon mass (Q 2 ~ .OlGeV 2 ) the integral is dominated by the low energy 
contribution to U R (Q 2 ). By doing the Taylor expansion of U R (Q 2 ) at Q 2 = and 
by imposing 11^(0) = one has 
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n£(Q 2 ) = qX'(Q 2 ) = q 2 

Notice that 11^' (Q 2 ) coincides with the first derivative only at Q 2 = 0. The term 
with the second derivative contains one additive power of the ratios Q 2 /A 2 and 
Q 2 /Mq, where A x ~ 1 GeV is the ultraviolet low-energy cutoff and Mq ~ 300 MeV 
is the infrared low-energy cutoff and are the natural dimensionful parameters of the 
long-distance expansion. 

By first approximation Ur(Q 2 ) ~ Q 2 ^r(0) and the integral of eq. (^) gives the 
value for afi, 

The calculation of reduces then to the calculation of the first derivative of 
the vector two-point function at Q 2 = 0, as was already pointed out in ref. @. In 
ref. [0 the long-distance behaviour of the vector function was derived in the ENJL 
framework. In |1J we have shown that the inclusion in the bosonized NJL Lagrangian 
of higher dimensional quark-resonance vertices, which are suppressed respect to the 
leading lowest dimensional four-quark operator by powers of the inverse cutoff A x , 
generates next-to-leading power corrections to the leading logarithms (NPLL) of the 
parameters of the effective resonance Lagrangian which are responsible of their Q 2 
dependence. 

As a consequence, higher dimensional quark-resonance operators modify the long 
distance behaviour of the vector Green's function predicted by the ENJL model 
and enter in the determination of through eq. (Q). Already in the first ap- 
proximation of eq. (0) the NPLL corrections proportional to Q 2 , i.e. of the type 
Q 2 /A 2 ln(A 2 /MQ) give contribution to the derivative at Q 2 = 0; the derivative at 
Q 2 = is sensitive to the behaviour at shorter distances. 

We proceed as follows. Using the first approximation of eq. (0) we review the 
ENJL prediction already derived in 0. Then we study the relevance of two sources 
of corrections: the NPLL contributions and the gluon contributions. Beyond the 
first approximation we analyze the sensitivity to both of them via the evaluation of 
the dispersive integral of eq. (f|) over the long-distance part of II^(Q 2 ). 

The vector two-point function in the ENJL model and in the chiral limit can 
be parametrized in terms of the running photon- vector coupling fv{Q 2 ) and the 
squared mass of the vector resonance My(Q 2 ) as follows: 



^R(O) + in 2 R (0) + 



(6) 
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UviQ ] ~ M 2 (Q 2 ) + Q 2 ' (8) 

or equivalently in terms of the vector two-point function in the mean field ap- 
proximation Uy(Q 2 ) [0 

c X 

where Gy(A x ) is the coefficient of the four-quark vector-like interaction in the 
ENJL Lagrangian. The Q 2 dependent parameters are given by: 



U v(Q 2 ) = C daa{l-a) T(0,a Q ) 

ion Jo 



fv(Q 2 ) = ^v(Q 2 ) 

with the incomplete Gamma function r(0, qq) = — lnctg — 7^ + 0(ocq), oq = 
(Mq + a(l — a)Q 2 )/A 2 and 7^ = 0.5772... is the Euler's constant. Expressions (§) 
and (||) with the parameters (|T0|) correspond to the diagram of Figure 2a, which 
is the infinite resummation of linear chains of constituent quark bubbles with the 
insertion of the leading four-quark vector operator of the ENJL Lagrangian. Linear 
chains of quark bubbles are of order N c , while loops of chains of quark bubbles are 
of order 1 in the 1/N C expansion. 

The derivative at Q 2 = is given by 

dQ 2 W " dQ 2 W [lvW N C A 2 X ~dQ* W M 2 (0Y { ) 

The ENJL model gives the following prediction: 



dTI* N c 1 4 , 51-<M Wn M ( 



dQ 2 '^ 16tt 2 M 2 15 



2 -, 

Q 



(12) 
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where (7,4 is the mixing parameter between the axial-vector and the pseudoscalar 
mesons in the bosonized ENJL action and it is related to the vector meson mass 
and the vector coupling Gy by the following relations [j/J: 
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iz^. 4 ^|r ( o,Mj/A^^. ( 13 ) 

With the values of the best fit 1 of ref. [§ M Q = 265MeV r , A x = 1.165GeV 
and qa = 0.61 we obtain for the value a!^ = 8.66 ■ 10~ 8 , which corresponds to the 
value ^t(O) = —0.164 of the first derivative. In formula ( |12"1) the incomplete Gamma 
function has been approximated to its leading logarithmically divergent part. This 
corresponds to the substitutions T(0,aQ) = — lnctQ — 7^ and exp(— Mq/A^) = 1, 
where exp(— Mq/A 2 ,) comes from the derivative of T(0, «q) at Q 2 = 0. 

Chiral loop corrections due to n, K exchanges, are not included in the above 
value of a^. They are given by the diagram in Figure 2b and are next-to- leading 
(i.e. 0(1)) in the 1/N C expansion; in the ENJL framework they are generated by 
diagrams with loops of chains of quark bubbles. This contribution has been derived 
in using ChPt with a value of (0.71 ±0.07) ■ 10~ 8 . The two summed contributions 
give 

a^had) = a^(Fig. 2a) + a M (xloops) = 9.37 • 10 -8 , (14) 

which is a rather high value compared to the phenomenological estimates in (|3|). 
Beyond the fact that the first derivative approximation can be not sufficiently accu- 
rate, we first analyze the effects of the two "extra" contributions we have mentioned. 

The first source of "extra" corrections can be derived in the framework of the 
Quark- Resonance model formulated in 

The vector two-point function calculated with the inclusion of NTL vertices cor- 
responds to the diagram of Figure 2c: the vector resonance exchange diagram plus 
the "local" diagram. The infinite resummation of linear chains of quark bubbles 
of Figure 2a is a part of the contributions to the renormalized vector resonance 
propagator, the insertion vector four-quark 

The complete set of NTL corrections (i.e. 1/A£), includes two types of contribu- 
tions U|: 

I. NTL genuine power corrections (NTLP) of the type: 

II. NTL power corrections to the leading logs (NPLL) of the type: 

o 2 i M o 1 

Class I is generated by an infinite number of higher dimensional quark-resonance 
vertices, while class II is generated by a finite number of 1/A^ vertices ||T|. A 
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best fit to the experimental e + e~ — ► hadrons cross section in the I = 1, J = 1 
channel and in the intermediate Q 2 region (500 MeV - 900 MeV) has shown the 
numerical relevance of the new NPLL counterterms proportional to Q 2 [p]], while 
NPLL corrections proportional to the IR cutoff Mq are negligible in this regime. 
They are not negligible when treating the very low energy behaviour (Q 2 — > 0) of 
the Green's functions, as it happens in the present case. 

In table (|IJ) the type of quark-resonance vertices we need are listed. We distinguish 
four sectors: the derivative one, the vector one, the scalar one and the scalar- 
vector one. NTL corrections proportional to Q 2 get contributions from the derivative 
and the vector sets, while NTL corrections proportional to Mq get contributions 
also from the scalar and scalar-vector sets when the scalar field assumes its VEV, 
< H >= Mq, which plays the role of the IR cutoff of the effective theory. 

The full running of the vector resonance parameters (Zy = vector wave function 
renormalization constant, fy = vector-photon coupling and My = vector mass) is 
parametrized as follows: 



z v (Q 2 ) 

M 2 {Q 2 ) 

MQ 2 ) 



N c 1 r 1 

— ^r- / da T(0, an) 



O 2 M 2 
6a(l - a) + I2a{\ - a)(3 v ^ + ^ M -^-Pi{a) 



A 2 

x 



1 



16tt 2 2G V ^ 
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Ml 
A 2 X Jo 



da T(0, aQ)P 2 (a) 



da T(0, an) 



} iQ 2 



6a(l - a) + 6a(l - a)/3^ + 

A x 



,i Q 2 



6a(l - a)(3 v — n , }l 



A 2 

x 



(15) 



In the calculation of fy and Zy we have not included NTLP corrections. They 
give a correction which is almost 1% of the leading logarithmic term. Genuine power 
corrections can be generated by the quadratically or more divergent part of diagrams 
that contain higher dimensional vertices which are suppressed by inverse powers of 
the UV cutoff A x . 

NTLP corrections to the squared vector mass and proportional to Mq can be 
reabsorbed in the renormalization of the vector coupling Gy as follows: 



A 2 

x 



G 



v 



A 2 

x 

Gy 



1 + 5' 



(16) 
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In eq. (|T5| ) the dependence upon the Feynman parameter a is always of the form 
a(l — a) for the Q 2 /A^ corrections. The polynomials -Pj(a) are explicitely calculable 
for each NPLL counterterm. We will put them equal to the leading polynomial 
a(l — a) which is still a good approximation. With this assumption the three vector 
parameters Zy, M v and f v at Q 2 = are given by: 



Z v (0) 
M 2 (0) 

fv(0) 



T&~r(o,M|/A- 



PI Ml 



A 2 

x 



16tt 2 2G V Z V (Q) 



9Z V (0) V16vr 2 
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6 A 2 



A x 



r(0,M 2 /A 
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(17) 



The NPLL corrections proportional to Mq do modify the ENJL leading prediction 
of the parameters of the effective meson Lagrangian at zero energy. One way to 
estimate the new coefficients (3 l M: is to do the best fit of the whole set of the leading 
(Mq, A x , qa) and NPLL parameters (with and without gluonic corrections) at Q 2 = 
using as inputs the experimental values of the low energy meson parameters; 
although a sufficiently accurate determination of the j3 l M is still not accessible with 
the present uncertainty on the very low energy experimental data. 

The ratio fy/My and the derivative of the squared coupling at Q 2 = 0, which 
enter in eq. (|TTD, are given by the following expressions and retaining up to 1/A 2 
terms: 



2 N r 



im 

M 2 (0) 316tt 2 M 2 (0) 



ENJL 



r(0,M 2 /A 



1 + 



P" m m) 



p 



A 2 

x 
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Again the approximation T(0, q;q) = — InctQ — ^ E and exp(— Mq/A 2 ) = 1 is 



ENJL 



understood and we have used the relation of eq. (|13|) which defines M v (0) 
assume its numerical value given by the fit 1 of ref. || which is (0.811) 2 GeV 2 . 
The size of the Q 2 corrections has been determined by the best fit of ref. ffl: 



We 



Pi 



-0.75 ±0.01 



p v = -0.79 ±0.01. 



(19) 
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The sign of /3p, (3 V is negative and increases the value of a^. By including only Q 2 
type corrections the value of increases to 1.22 ■ 10~ 7 . Comparable values of the 
(3 l M coefficients give as a maximum range of variation of aj 1.20 • 1(T 7 + 1.24 • 10~ 7 . 
This proves that Q 2 dependent contributions give the bulk of the NPLL corrections 
to calculated within the first derivative approximation. 

Gluonic corrections can be parametrized following ref. 0. The leading contribu- 
tion in the 1/N C expansion involves only one unknown parameter g which is related 
to the lowest dimensional gluon vacuum condensate: 



n 2 



x , , M < —GG > . (20) 

Because < ^-GG > is 0(N C ) g is 0(1) and the leading gluonic correction to 
the two-point vector function is still 0(N C ). The role of gluonic corrections in 
effective low energy fermion models is still an unsolved theoretical problem. What 
is clear is that the two gluon condensate of eq. (|20"D is only the low energy (< 
A x ) 'residue' of the standard two gluon condensate which is phenomenologically 
estimated through QCD sum rules. The latter suggest an 0(1) g parameter, while 
best fits of the ENJL parameters, using as inputs the experimental values of fir,Li 
and the meson resonances' parameters ||, strongly favour a value of g < 0.5. The 
first phenomenological estimation, usually referred to as the "standard value" , has 
been obtained by Shifman et al. (SVZ) J| by studying the charmonium channel 
and using Operator Product Expansion: they obtain < ^-GG >= 0.012 GeV 4 (i.e. 
g=1.3). A compatible value < a s GG >= (3.9 ± 1.0)10~ 2 GeV 4 has been obtained 
||10"| from the e + e~ — > / = 1 hadron cross section and using moment sum rules 
ratio. The most recent estimation via FESR (Finite Energy sum rules) JD] gives a 
significative higher value < ^GG >= 0.044±g;g^ GeV 4 (i.e. g ~ 4.8), although the 
error in all cases has to be conservatively taken around 40%. 

In what follows we will use three values of g which are acceptable in the ENJL 
model. The IR cutoff Mq, the UV cutoff A x and the axial-pseudoscalar mixing 
parameter g^, extracted from the best fits of experimental data, vary as functions 
of g as it is summarized in table (||). Mq decreases sensitively by increasing g. 

Leading gluonic corrections can be expressed as an additive contribution to the 
function U V (Q 2 ) which is written for N c = 3 as |12] 



n^(Q 2 ) = |^(-l + 3J 2 -2J 3 ), (21) 
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in terms of the Q 2 dependent functions 



f 1 1 
In = da m . (22) 



Expanding at small Q 2 we obtain their values at Q 2 = 0: 



n^(o) = -- ^ n'^(o) = — (23) 

y v ; 20vr 2 v v ; 70n 2 M? ' 



which give the following expression for the derivative: 



^V(O) = ^ (Q) + A g 
dQ 2K) dQ 2 y) 



69 PM 2 ^^^ 



70tt 2 M£ V20tt 2 / iV c A 2 V 3# 



1 „_Q J <? — 

Notice that Uy (0) = II ^ (0); all the quantities with superscript g = are 
those at g=0 with the parameters Mq, A x and rescaled according to table (§) 
for a given g in formula (23). In table (|3|) we give the gluonic corrections to the 



derivative at Q 2 = with g = 0.25, 0.5. 

Gluonic contributions decrease towards a better agreement with the phe- 
nomenological estimates ([|). 

The evaluation of the full dispersive integral (f|) requires the knowledge of the 
long distance (Id) plus the short distance (sd) behaviour of Il^Q 2 ). We can do our 
best performing the matching between the Id prediction coming from the effective 
theory and the sd prediction coming from perturbative QCD. In ref. [0] a value of 
a^(Fig.2a) = 6.7- 10~ 8 has been obtained in the ENJL framework with a best fitted 



matching point x ~ 0.91 which corresponds to an euclidean Q = x /(l — x)m 
(320) 2 MeV 2 . The value obtained for is quite better than the first approximation 
value 8.66 • 10" 8 . 

To see how "extra" corrections modify the ENJL prediction is sufficient to study 
the integral (|j) over the Id part in the range Jq for different values of x corresponding 
to an equivalent value of Q 2 = (0.3) 2 , (0.5) 2 , (0.8) 2 GeV 2 . 

NPLL corrections proportional to Q 2 in the QR model lead to a vector two-point 
function which decreases faster in Q 2 then the ENJL prediction. The dispersive 
integral gives for the values of table (|j). 
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The correction induced by higher order terms proportional to Q 2 and which are 
relevant in the intermediate Q 2 region is about one percent. This proves that is 
practically only sensitive to perturbative corrections which modify the very low Q 2 
region, i.e. Q 2 < (500MeV) 2 . 

Gluonic corrections modify the ENJL vector two-point function for all Q 2 as 
shown in figure 3. They modify the ENJL prediction of as summarized in table 
(|J). The corrections are 10% for g=0.5 and 15% for g=0.25. They decrease all. A 
better determination of non-gluonic contributions to can be used to constrain the 
value of the g parameter in low energy effective fermion models. 

We conclude that the ENJL prediction is reliable within 30%. Both gluonic cor- 
rections and next-to-leading higher dimensional quark-resonance interactions have 
to be taken into account if one wants to reach an accuracy better than 30%. The 
first derivative approximation is not sufficiently accurate to estimate the hadronic 
vacuum polarization contribution to a^, although the quantity is sensitive only to 
"extra" corrections in the low Q 2 region (Q < 500 MeV). NPLL corrections propor- 
tional to Q 2 do not affect a^. NPLL corrections proportional to Mq and gluonic 
corrections are relevant in the low Q 2 region. Their inclusion can explain the full 
agreement of the ENJL prediction with the phenomenological estimates @ of a 1 ^. 



Acknowledgements. I am grateful to Eduardo de Rafael for having called my 
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FIGURE CAPTIONS 

1) Hadronic vacuum polarization diagram to the anomalous magnetic moment of 
the muon. 

2a) The vector two-point function in the ENJL model and leading in the 1/N C ex- 
pansion. It is given by the infinite resummation of linear chains of constituent 
quark bubbles; each four-quark vertex is the leading vertex of the ENJL model 
with coupling Gy. 

2b) Chiral loops corrections to the vector two-point function. They are of 0(1) 
(i.e. next-to-leading) in the 1/N C expansion and in the ENJL model they are 
generated by the infinite resummation of loops of chains of constituent quark 
bubbles. 

2c) The vector two-point function in the QR model and leading in the 1/N C expan- 
sion. It is given by the "local" diagram and the vector-exchange diagram with 
the renormalized vector meson propagator. The diagram of figure 2a is a part 
of the contribution to the renormalized vector meson propagator. 

3) The vector two-point function U V (Q 2 ) in the ENJL model without gluonic cor- 
rections (solid line) and with gluonic corrections for g=0.25, 0.5 (dashed lines). 



TABLE CAPTIONS 

1) 1/A^ quark-resonance vertices of the sectors: I) derivative, II) vector, III) scalar 

and IV) scalar-vector which give contribution to the NPLL corrections. 

2) Values of Mq, A x and g& of the ENJL model obtained from the best fits of ref. 

|| to the experimental data of low energy parameters and with fixed values 
of the gluonic paramater g=0, 0.25, 0.5. 

3) Gluonic corrections to cfc in the first derivative approximation for two values 

of the gluonic parameter g favoured by the ENJL model. A 9 is the gluonic 
correction to the first derivative as defined in eq. (|24j), dll v / dQ 2 (0) is the first 
derivative including gluonic corrections. In the last two columns the numerical 
value of a! 1 and the variation in percentage are shown. 
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4) Numerical values of obtained through the dispersion relation (^), where the 

integral is performed on the long-distance part of IT^, predicted by the QR 
model in the range < Q 2 < Q 2 and compared to the ENJL prediction. 

5) Numerical values of obtained through the dispersion relation (^), where the 

integral is performed on the long-distance part of IT^. predicted by the ENJL 
model in the range < Q 2 < Q 2 without the inclusion of gluonic corrections 
(g=0) and including gluonic corrections with g=0.25, 0.5. 
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DERIVATIVE 

ft QludPT^Q + ft Q^{7, t}Q 
VECTOR 

ftQl^W+Q + ft v Qi,{W+, t}Q + ftQ^{(d„ Zh W+}Q + ftQi,[V^ W+}Q+ 
ftQl,{W +2 , 7}Q + ftQi^Wf, W+]Q + ftOy^W+W? d„ + Z W+W+)Q+ 
ftQl,{W+W+ Z + Z W+W+)Q + ftQ-y^W;, W+]Q 

SCALAR 

ft QH 3 Q + ft Q^{H 2 , 7}Q + ft Q{H> ~d}Q 
SCALAR-VECTOR 

ft v Q{W +2 , H}Q + ft v QW+HW+Q + ft v Q^{H 2 , W+}Q + ft v Q llx HW+HQ 

+ft v Q(w;H d, + X HW Z)Q + PIvQ(hw+ d, + d, w;h)q. 

Table 1: 

The covariant derivative of the vector field Wij" and the scalar field H is defined 
as d^O = dpO + [r M ,C?], while the covariant derivative on the constituent quark 
field Q is defined as d^Q = d^Q + T M Q. = 1/2(^8^ + £8^) is the vector 
current constructed with the square root of the pseudoscalar meson field U = £ 2 = 
exp(i(f) I f n ) . The identity = —\f£ v + introduces the field strenght of 

the electromagnetic field ft, = £F^ U ^ + which generates the photon-vector 

interaction associated with the coupling f v . 
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Mq (MeV) 


A x (GeV) 


9A 


g = o 


265 


1.165 


0.61 


g = 0.25 


246 


1.062 


0.62 


4 = 0.5 


204 


1.090 


0.66 



Table 2: 



9 


A 9 


dU l v /dQ 2 {0) 


< 


Aa,/a,(%) 


g = 0.25 


0.048 


-0.136 


7.2 • 10~ 8 


-17% 


9 = 0.5 


0.134 


-0.133 


7.0 • 10~ 8 


-19% 



Table 3: 



Q 2 GeV 2 


a® R x 10 s 


a EJVJL x 1Q 8 


(0.5) 2 
(0.8) 2 


6.9 
7.3 


6.8 
7.1 



Table 4: 
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Q 2 GeV 2 


<*5 


g=0 




=0.25 




;=0.5 


(0.3) 2 


5.9- 


io- 8 


5.0 


■ io- 8 


5.3 


io- 8 


(0.5) 2 


6.8- 


io- 8 


5.9 


■ io- 8 


6.3 


io- 8 


(0.8) 2 


7.1 ■ 


10~ 8 


6.2 


■ io- 8 


6.6 


io- 8 



Table 5: 
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